co = Crystal.create_hematite_rhombohedral;
%co.convert_rhombohedral_to_hexagonal;

clf;
co.plot_atoms;

pos = POSCAR('hematite', 1, co.LatticeVectors,...
        co.NumAtoms, 'C', co.Positions);
pos.write;